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Studies of rare B hadron decays to leptons at hadron colliders 
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Abstract. Rare B hadron decays provide an excellent test bench for the Standard Model and can probe new 
physics models. We review the experimental progress of the searches for rare leptonic B decays (b — > and 
b — > sf + £~) at LHC and Tevatron experiments. 



1 Introduction 

B hadron decays mediated by flavor-changing neutral cur- 
rents (FCNC) are suppressed in the Standard Model (SM) 
and are sensitive to new particle contributions through 
loop diagrams. Searches for B" d meson decays to dimuons 
were performed both at the Tevatron and the LHC, where 
the large samples of B hadrons provide a unique opportu- 
nity to study rare decays. In addition, the rates and angular 
distributions of b — > sp + p" decay products are also sensi- 
tive to new physics. In Section [2] we review the experi- 
mental status of the searches for the B° and B° rare decays 
to dimuons at hadron colliders. In Section [3] the angular 
analysis and measurements of branching ratios for various 
b — > sp + p~ decays performed by the CDF experiment is 
summarised. 

2 Searches for the rare decays B° d -> pi^pc 

The rare decays B^ — > p + p and B° — > p + p~ are highly 
suppressed in the SM of particle physics due to their 
FCNC nature. These decays are forbidden at tree level 
and can proceed only through higher order diagrams, such 
as electroweak penguin and box diagrams. Furthermore, 
the decays are helicity suppressed and require an internal 
quark annihilation within the B meson. The SM predic- 
tions for the B° and B° branching ratios to dimuons are 
(3.23 + 0.27) x and (1.07 + 0.10)x lfr 10 , respectively, 
with the main uncertainty resulting from the value of the B 
meson decay constant f B [ 1 ]. The comparison with exper- 
imental results requires the inclusion of soft-photon radia- 
tion and B" -B° s oscillations effects each yielding 0(10%) 
corrections to the predicted decay rate. Several extensions 
of the SM, such as supersymmetric models and models 
with a non-standard Higgs sector, predict enhancements 
or suppressions to the branching fractions for these rare 
decays. The rare nature of the processes and the rather pre- 
cise SM predictions make these decays an excellent probe 
for physics beyond the SM. 
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At the LHC searches for the Bj, and B decays to 
dimuons were performed by the ATLAS, CMS and LHCb 
experiments using datasets of pp collisions at center-of- 
mass energies -\fs — 7 or 8 TeV collected between 2010 
and 2012. The CMS experiment performed a simultane- 
ous search for the rare decays B'' — > p + p~ and B° — > /lo- 
using an integrated luminosity of 5 ftr 1 collected at the 
center-of-mass energy yfs -1 TeV @. An event-counting 
experiment was performed in dimuon mass regions around 
the B ( s * and B° masses and all selection criteria were estab- 
lished before observing the signal region. A normaliza- 
tion sample of events with B + — > J/i//K + decays (where 
J/t// — > p + p~) was used to remove uncertainties related 
to the b-quark production cross section and the integrated 
luminosity. Combinatorial backgrounds were evaluated 
from the data in dimuon invariant mass sidebands while 
backgrounds from B decays were assessed with Monte 
Carlo (MC) simulation. The analysis was performed sep- 
arately in two channels, barrel and endcap, and then com- 
bined for the final result. The barrel channel contained 
the candidates where both muons have I77I < 1.4 and the 
endcap channel included those where at least one muon 
had I77I > 1.4. B— > yu + yu~candidates with traverse mo- 
mentum above 6.5(8.5) GeV in the barrel(endcap) were 
formed by two oppositely-charged muons originating from 
a common vertex and with an invariant mass in the range 
4.9 < m Mfl < 5.9 GeV. The (sub)leading muon transverse 
momentum was required to be larger than (4.0)4.5 GeV 
in the barrel and (4.2)4.5 GeV in the endcap. Further 
selection cuts were applied on the B candidate isolation 
variables, decay length significance and three-dimensional 
pointing angle. 

The branching fraction was measured using the follow- 
ing equation 

S(B ^ p V) = — - — S(B + -> J/^K + ), (1) 

norm fs ^sig 

where e S i g (e, wrm ) is the total signal(normalization) effi- 
ciency, N, wrm is the number of reconstructed B + — > J/tfrK + 
decays, S(B + ) is the branching fraction for the normaliza- 
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Figure 1. Dimuon invariant-mass distributions measured by CMS in the barrel (left) and endcap (right) channels. The signal windows 
for B° and B° are indicated by horizontal lines 



Table 1. Expected and observed upper limits (95% CL) for the B s d branching fractions to dimuons. 
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tion channel, f u /f s is the ratio of the B + and B' produc- 
tion cross sections, and N S j g is the background-subtracted 
number of observed candidates in the signal window. Fig- 
ureQ~]shows the measured dimuon invariant-mass distribu- 
tions. Six events were observed in the B" — > // + /LTsignal 
windows, while two events were observed in the B° — > 
n channel. This observation is consistent with the SM 
expectation for signal plus background. Exclusion lim- 
its on the branching fractions were obtained with the CL S 
method and are reported in Table [T] 

The ATLAS experiment performed a search for the 
rare decay B'' — > /i + yu~using 2.4 fb -1 collected at V* = 
7 TeV Q. The sensitivity to the B° — > yU + |/~decay is be- 
yond the reach of the current analysis. Hence only a limit 
on the B° decay was derived by assuming the B° branching 
ratio to be negligible. The di-muon mass region 5.066 < 
trifift < 5.666 GeV was removed from the analysis until the 
procedures for event selection, signal and limit extractions 
were fully defined. The sample of signal candidates was 
selected with a multivariate classifier, trained on a frac- 
tion of the events from the di-muon invariant-mass side- 
bands. 14 discriminating variables were used in the selec- 
tion, including the isolation, decay angle and proper de- 
cay length significance, vertex separation, impact param- 
eter of the decay products as well as the B hadron and 
single muon transverse momenta. Three regions of differ- 
ent mass resolution and hence signal-to-background ratio 
were defined and selection cuts were optimized indepen- 



dently. The three categories were defined by the intervals 
Vlmax\ = 0—1, 1 — 1.5 and 1.5-2.5, where rj max is the largest 
pseudorapidity value of the two muons in each event. The 
width of the search region (116 to 171 MeV) and the clas- 
sifier output threshold (0.234 to 0.270) were optimised fol- 
lowing the prescriptions of Ref. [4|. The MC validation for 
the kinematic distributions and classifier output was per- 
formed on the normalization channel B + — > J/i[/K + and 
residual discrepancies between data and simulation were 
treated as systematic uncertainties. The independence of 
the classifier output on the dimuon mass was tested by re- 
training the same classifier in an higher-mass unblinded re- 
gion. The branching ratio was obtained with Eq. [1] where 
the acceptance and selection efficiencies were extracted 
from MC simulation. In each mass-resolution category 
the Bj — » ix+iT signal yield was obtained from the number 
of events observed in the search window (see Figure 131 , 
the number of background events in the side-bands (ex- 
cluding the events used in the classifier training and cut 
optimisation processes), and the small amount of resonant 
background. The expected and observed upper limits on 
the branching ratio at 95% confidence level (CL) were ob- 
tained with the CL S method and are reported in Table [1] 

The exclusion limits on the branching ratios for 
B'' and B° decays to dimuons obtained from the ATLAS, 
CMS and LHCb searches were combined using the CL S 
method 0. The LHCb results used in the combination 
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Figure 2. Invariant-mass distribution of dimuon candidates 
in ATLAS data for the pseudorapidity region with the best 
invariant-mass resolution. The invariant-mass distribution for the 
selected candidates in data (dots), the signal (continuous line) as 
predicted by MC assuming S(B° — > yu + yu")=3.5 x 10" 8 , and two 
dashed vertical lines corresponding to the optimized signal re- 
gion. The grey areas correspond to the sidebands used in the 
analysis. 



were obtained analyzing integrated luminosities of 1 fb 1 
and 0.037 fb" 1 collected in 201 1 and 2010, respectively, at 
yfs = 7 TeV iQ. No ATLAS data were used in the 
limit combination for the B° decay. The results of the 
LHC combination are reported in Table [T] and constitute 
the most stringent exclusions limits to date for both de- 
cays. 

The LHCb experiment recently updated the search 
combining datasets corresponding to integrated luminosi- 
ties of 1.1 fb -1 and 1.0 fb -1 collected at yfs = 8 TeV and 
7 TeV, respectively |8j. The analysis yielded an excess 
of events in the — > yu + //"channel with respect to the 
background expectation, with a signal significance of 3.5 
standard deviations. An unbinned maximum-likelihood fit 
gave a branching fraction S(B° -> p + pT)= (3.2+}|)x10" 9 , 
where the statistical uncertainty is 95% of the total uncer- 
tainty. This value is in agreement with the SM expecta- 
tion. The number of candidates in the B° — > /^^"mass 
window is consistent with the background expectation 
and the observed(expected) 95% CL exclusion limit is 
9.4(7.1) x 10" 10 . A more detailed discussion of the LHCb 
results can be found in these proceedings |9 1. 

At the Tevatron similar searches were performed by 
the CDF and DO experiments. The most stringent con- 
straint is provided by the DO collaboration, that has fully 
analyzed the Run II pp dataset corresponding to an inte- 
grated luminosity of 10.4 ftr 1 IfTUl . The mass resolution 
is not sufficient to separate the B" and B° decays and a 
search for decays was performed assuming no contribu- 
tion from B° . Decay candidates were identified by select- 
ing two high-quality muons of opposite charge forming a 
three-dimensional vertex well separated from the primary 
interaction. Backgrounds in the low mass region are dom- 
inated by sequential decays b — > cpX —> yiyfii' while on 
the high mass region double semileptonic decays are dom- 
inant. From the peaking backgrounds the B° — >KK is the 
dominant contribution. Oppositely charged muons with 



Pt > 1.5 GeV were selected and a multivariate discrim- 
inant was used to differentiate between signal and back- 
grounds. Two boosted decision trees (BDT) were trained 
against the two background categories using 30 variables. 
The optimal BDT output cuts were determined by opti- 
mizing the expected limit the B 1 ' branching ratio and were 
set to 0. 19 and 0.26. The expected number of events in the 
signal and control regions was determined by applying a 
log likelihood fit to the dimuon mass distribution using an 
exponential plus constant functional form. Three events in 
the signal region were observed after unblinding, which is 
consistent with expected background. The observed and 
expected upper limits improve substantially the previous 
DO result and are given in Table [T] 

The CDF collaboration has analyzed the full Run II 
data sample ifTTl . with an event selection based on a neu- 
ral network using 14 input variables. It required high 
quality muon candidates with transverse momentum p T > 
2.0(2.2) GeV in the central (forward) region. The muon 
pairs were required to have an invariant mass in the range 
4.669 < m Ml j < 5.969 GeV and are constrained to originate 
from a common well measured three-dimensional vertex. 
A likelihood method together with a selection based on the 
track energy loss were used to further suppress contribu- 
tions from hadrons misidentified as muons. The combi- 
natoric background was estimated by fitting a fixed slope 
first order polynomial to the mass sidebands while the 
peaking background from B— > hti decays is evaluated 
from MC, where the hadron misidentification probabili- 
ties were extracted from data. The number of candidates 
in the B° mass window after unblinding are consistent 
with the background expectations and the corresponding 
exclusions limits are reported in Table [T] The search in 
the B" window yielded an excess of events in the region 
of neural network output v > 0.97 with a background- 
only p-value of 0.94%. The combination of all neural net- 
work bins has a p-value of 6.8% with respect to a back- 
ground plus SM-signal hypothesis. The corresponding 
observed and expected exclusion limits obtained with the 
CL S method are given in Table [T] 

The impact of the B s d— > //^"searches at hadron 
colliders on new physics models is represented in Fig- 
ure]^ showing the correlation between B° and B° branch- 
ing ratio to dimuons for models with Minimal Flavor 
Violation (MFV), the SM with a sequential fourth gen- 
eration (SM4) and four SUSY flavour models (MSSM- 
RVV2/AKM/AC/LL) (T2|. A large part of the parameter 
space of the supersymmetric models, where tan (3 can be 
large, is ruled out by the current constraints. While the re- 
cent LHCb observation indicates a value of S(B" — > [i + [T) 
consistent with the SM, large enhancements of S(B° — > 
p + p~) are not yet excluded by experimental data. The most 
stringent exclusion limit to date for this channel is about 
eight times larger than the SM prediction (at 95% CL). The 
prospects for further improving the precision on the B ! £ / 
branching fractions to dimuons at the LHC are excellent. 
The LHCb, CMS and ATLAS experiments have collected 
about 2.1 fb" 1 , 23 fb" 1 and 22 fb" 1 in 2012, respectively, at 
the center-of-mass energy yfs = 8 TeV. Results based on 



EPJ Web of Conferences 



the analysis of these datasets are expected to be released 
within 2013. 
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Figure 3. Summary of the measured exclusion limits at 95% 
CL for the B° and B° decays to dimuons at LHC and Tevatron. 
Areas on the right hand side of the dashed line are excluded. The 
regions predicted by SM4 and four SUSY models are represented 
by the colored contours. The SM expectation is marked by a star 
(from fl2l ). 
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Rare B hadron decays of the type b — > s/i + /j. are medi- 
ated by FCNC and occur in the SM with branching ratios 
of the order O(10~ 6 ) [13] [14). Theories extending the SM 
favour enhanced rates of these decays and can modify the 
angular distributions of the decay products. Although no 
significant departures from the SM have been observed so 
far, the LHCb collaboration recently reported an isospin 
asymmetry between the branching ratio of B + — > K + /i + jiT 
and B° — > K°^ + /U deviating from zero with a significance 
of 4.4 standard deviations lfT51 . In addition, the Belle ex- 
periment has measured a lepton forward-backward asym- 
metry (A FB ) different from SM expectation at 2.7 standard 
deviations 1 16 1. 

The CDF experiment recently analysed the full dataset 
corresponding to an integrated luminosity of 9.6 fbr 1 . Be- 
sides the measurements of total and differential branch- 
ing ratios in B + -> K + jU + j"~, B° ~> K*(892)VV ; 
B° -> <^< V . B° -> KVV , B + -> K*(892) + yuV, and 
A° — > Kn + n~, measurements of the combined branching 
ratio assuming isospin symmetry and of the isospin asym- 
metry between neutral and charged B mesons are pro- 
vided [17|. To cancel the dominant systematic uncertain- 
ties, the decay rates for each rare channel Hb — > hyfpT 
were normalised to the corresponding resonant channel 
H b -> hJ/iff, where H b can be B + , B°, B° and A b , and 
h stands for K + , K*°, 0, K°, K* + , and A. Candidate events 
were selected by constructing a vertex of two muons that 
satisfy the trigger requirements with one charged track or 
with two reconstructed tracks of opposite charge, for the 
case of charged and neutral h hadrons, respectively. For 
the normalisation samples the dimuon invariant mass was 
required to be within 50 MeV of the J /if/ mass. After loose 
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Figure 4. Up: Invariant mass distribution for the B° — > K*°;u + /r 
channel measured by the CDF experiment. Down: Differential 
branching ratio as function of q 2 (squares) compared to the SM 
expectation (red line). 



selection cuts the rare decays were tightly selected with 
a multivariate analysis based on a Neural Network. The 
signal yields were obtained from an unbinned maximum 
likelihood fit to the invariant mass distribution, as shown in 
Figure |4] and were corrected for the MC selection efficien- 
cies to extract the relative branching ratios. The results are 
summarised in Table [2] The differential branching ratios 
were measured as function of the dimuon invariant mass, 
q 2 = m^c 2 , by performing the signal fit in each q 2 bin, as 
shown in Figure 0for the B° -> K* W decay. 

The isospin asymmetries between neutral and charged 
B meson decays are observables with relatively low theo- 
retical uncertainties, since the leading form factor uncer- 
tainties cancel in the ratio. The asymmetries 



S(B° -> K W VV) -S(B + -> K w V/0 
S(B° -> K^Vl^) + S(B+ -> KW + /i+/i-) 



were also measured by the CDF experiment and were 
found to be consistent with zero over the full q 2 range. 
Thus this result does not confirm the deviation observed 
by LHCb. 

In B° — > K ( * ) /i + ; u _ decays the muon forward-backward 
asymmetry (Afb) and the K* longitudinal polarisation (Fl) 
are extracted from the angular distributions cosfl^ and 
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Figure 5. Longitudinal polarisation (left) and forward-backward asymmetry (right) for B° — > K*/u + fi decays as function of the dimuon 
invariant mass. The SM expectation is represented by the solid line. 



Table 2. Relative branching ratios for b — > s/j + fi decays measured with the CDF experiment. 



Channel ratio 


Relative branching ratio (10 3 ) 


S(B + - 


. KV/0/£(B + 


— » 


K + 7/-A) 


0.44 ± 0.03(stat)±0.02(syst) 


S(B° -> 


K*VV)/S(B° 


— > 


K*°J/l//) 


0.85 ± 0.07(stat)+0.03(syst) 


S(B° 




— > 




0.90 ± 0.14(stat)±0.07(syst) 


S(B° - 


> kVV)/^(b° 


— > 


K°J/tff) 


0.38 ±0.10(stat)±0.03(syst) 


S(B + — > 


K* + /i + fT)/8(B + 


— > 


K* + J/ifr) 


0.62 ± 0.18(stat)+0.06(syst) 


£(A°- 


-» A^V)/B(A° 




AJ/i/f) 


2.75 + 0.48(stat)+0.27(syst) 



cos 9 K , respectively, where 6^ is the angle between the 
direction and the opposite of the B meson direction in the 
dimuon rest frame, while 9k is the angle between the kaon 
direction and the direction opposite to the B meson in the 
K* rest frame. A simultaneous unbinned maximum likeli- 
hood fit of the polarisation angles is performed to extract 
the angular variables as function of the dimuon invariant 
mass, as shown in Figure [5] All measurements are con- 
sistent with the SM prediction and with previous measure- 
ments. 
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